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ABSTRACT 

A  horizontal  circular  cylinder,  located  near  the  floor  of  a  two 
dimensional  wave  channel,  was  subjected  to  a  train  of  gravity  waves.  The 
wave  height  and  horizontal  and  vertical  forces  on  the  cylinder  due  to  the 
incident  waves  were  measured  and  presented  in  dimensionless  form.  Experi- 
mental values  of  the  horizontal  and  vertical  dimensionless  force  coeffi- 
cients are  presented  as  functions  of  relative  wave  height  and  relative 
wave  length.  The  dimensionless  force  coefficients  predicted  by  a  modified 
Morrison's  equation  are  compared  with  the  experimental  data.  Experimental 
wave  lengths  of  from  two  to  seventeen  feet  were  investigated. 

The  horizontal  force  coefficients  were  found  to  vary  linearly  with 
wave  height  over  the  range  of  wave  lengths  tested.  The  vertical  force 
coefficients  displayed  regions  in  which  the  force  was  both  inertia 
dominated  and  dominated  by  a  lift  force. 
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I.   INTRODUCTION 

Exploration  of  the  ocean  floor  has  received  increased  interest  in  the 
past  several  decades,  particularly  the  area  of  the  continental  shelf.  A 
significant  amount  of  attention  has  been  devoted  to  offshore  oil  explo- 
ration, to  the  recovery  of  petroleum  and  to  deployment  of  undersea 
habitats.  Often  petroleum  is  located  in  remote  areas  of  the  world  and 
oil  production  involves  use  of  very  large  submerged  oil  storage  tanks  and 
associated  pipelines.  In  other  locations  petroleum  has  been  piped  ashore, 
through  the  surf  zone,  to  shore  facilities.  Recently,  interest  has  been 
demonstrated  [1]  in  the  deployment  of  large  submerged  nuclear  power 
plants.  In  the  reference  cited  it  is  proposed  that  a  nuclear  plant  be 
placed  on  the  ocean  floor  in  two  hundred  and  fifty  feet  of  water. 

As  one  might  expect,  this  activity  has  generated  considerable  interest 
in  the  interaction  of  gravity  waves  with  submerged  objects.  For  example, 
the  wave  force  exerted  on  vertical  piling  has  been  a  topic  of  many  investi- 
gations during  the  past  twenty  years  [2]  ,  [3] ,  [4]  .  In  these  studies  the 
well  known  Morrison  Equation  [4]  appears  to  be  the  most  practical  approach 
to  the  calculation  of  wave  forces.  Ihis  equation  combines  both  a  drag  and 
inertia  term,  the  values, of  which  must  be  determined  experimentally.  At 
best,  these  two  experimental  coefficients  are  dependent  upon  the  geometry 
of  the  pile  and  amplitude  of  the  fluid  motion,  and  it  would  appear  that 
twenty  years  of  research  has  only  partially  delineated  this  interdependence, 

In  the  case  of  large,  deeply  submerged  objects  the  Morrison  Equation 
is  also  applicable  and  generally  yields  valid  results  for  the  horizontal 
wave  forces  provided  diffraction  effects  do  not  become  appreciable.  In 
this  instance,  the  drag  term  may  be  discarded  since  the  hydrcdynamic  force 
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is  primarily  inertial  in  nature.  However,  it  is  still  necessary  to  know 
the  added  mass  coefficient  as  it  depends  on  the  body  configuration  and 
the  effect  of  the  rigid  bottom  surface.  Moreover/  if  the  object  is 
bottom  mounted  the  vertical  fluid  acceleration  is  zero,  and  accordingly, 
the  Morrison  Equation  yields  zero  vertical  force,  a  result  that  is 
obviously  invalid. 

A  diffraction  theory,  valid  for  the  calculation  of  forces  exerted  on 
large  submerged  or  semi-submerged  bodies  of  arbitrary  shape,  which  accounts 
for  both  the  effect  of  the  free  surface  and  bottom,  has  been  developed  by 
Garrison  and  Chow  [5]  .  Theoretical  results  were  compared  with  corres- 
ponding experimental  results  for  two  different  practical  submerged  tank 
configurations  and  good  agreement  was  found.  However,  the  method  is  valid 
for  linear  incident  waves  only;  the  nonlinear  effects  of  large  amplitude 
waves  can  not  be  determined. 

In  the  case  of  large  objects,  which  are  not  so  deeply  submerged  that 
the-  free  surface  has  a  sizable  effect,  no  theory  is  available  which  is 
valid  for  large  values  of  wave  height  and  few  experimental  results  have 
been  published.  However,  one  study  has  been  carried  out  by  Johnson  [6] 
for  horizontal  wave  forces  acting  on  a  bottom-mounted,  horizontal 
circular  cylinder  in  long  waves.  Inasmuch  as  only  long  waves  were  con- 
sidered, the  results  were  dependent  on  the  wave  height  and  water  depth 
only  and  no  consideration  was  given  to  intermediate  wave  lengths.  More- 
over, rather  than  correlate  the  forces  with  appropriate  dimensionless 
parameters,  Johnson  used  a  regression  analysis  to  correlate  the  forces. 
Consequently,  little  light  was  shed  on  the  basic  mechanism  involved. 


Schiller  [7]  measured  wave  forces  acting  on  a  submerged  horizontal 
circular  cylinder  due  to  rather  small  amplitude  waves.  Over  the  range  of 
wave  heights  considered  in  his  tests,  he  found  the  amplitude  of  the  hori- 
zontal force  to  vary  linearly  with  the  wave  height.  At  small  wave  heights 
the  vertical  force  was  found  to  vary  harmonically  with  the  wave  motion, 
the  maximum  amplitude  of  which  increased  linearly  with  the  wave  height. 
At  larger  values  of  wave  height,  the  vertical  force  amplitude  tended  to 
increase  in  proportion  to  the  wave  height  squared.  This  latter  variation 
of  the  wave  force  is  caused  by  the  velocity  squared  term  occuring  in 
Bernoulli's  Equation. 

The  present  investigation  was  an  extension  of  Schiller's  work  to 
include  large  wave  heights.  However,  the  present  investigation  was  limi- 
ted to  the  single  configuration  where  the  horizontal  circular  cylinder 
was  placed  on  the  bottom.  A  transition  section  was  placed  in  the  wave 
channel  used  by  Schiller,  at  a  location  twenty  feet  from  the  wave 
generator  for  purposes  of  providing  a  natural  transition  from  intermediate 
or  deep  water  waves  to  shallow  water  waves  of  finite  amplitude.  The  beach 
slope  was  adjusted  for  negligible  reflection  as  discussed  in  the  section 
on  experimental  apparatus. 
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II.   THEORETICAL  ANALYSIS 

The  problem  under  consideration  is  depicted  in  Figure  1.  A  train  of 
regular  waves  of  height  H  is  considered  to  progress  in  the  positive  x 
direction  in  water  of  depth  h.  It  is  the  present  interest  to  determine 
the  horizontal  and  vertical  components  of  wave  force  acting  on  the 
horizontal  circular  cylinder  in  contact  with  the  rigid  bottom. 

The  exact  analytical  solution  to  this  problem  is  extremely  difficult 
and,  therefore,  only  an  approximate  analysis  is  attempted.  However,  it 
is  first  instructive  to  carry  out  a  dimensional  analysis  of  all  of  the 
pertinent  parameters.  It  is  known  a.  piiosii   that  the  maximum  horizontal 
or  vertical  force  per  unit  length  of  a  cylinder  is  dependent  upon  the 
following  variables: 

x  max    y  max   -,.   T  TT     N  ,,» 

or  -^ =  f  (h,a,L,H,p,g,y)  (1) 


in  which 


H  =  wave  height 

I   =  cylinder  length 

h  =  water  depth 

a  =  cylinder  radius 

L  =  wave  length 

p  =  fluid  density 

g  =  gravitational  acceleration 

y  =  fluid  viscosity. 


It  may  be  noted  at  this  point  that  a  relationship  exists  between  the 
parameters  associated  with  the  incident  wave  (i.e.,  H,  T,  L  and  h,  where 
T  denotes  the  period) .  Consequently,  only 'three  of  the  incident  wave 
parameters  are  needed  in  the  dimensional  analysis,  i.e.,  either  H,  h  and 
L  or  H,  h  and  T.  Moreover,  either  of  these  three  sets  of  parameters 
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FIGURE  1  DEFINITION  SKETCH 
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completely  describes  the  incident  wave  regardless  of  whether  or  not  it 
is  of  small  amplitude.  For  small  amplitude  waves  the  well-known 
relationship 

2JL  =lTanh2lh  (2) 

gT2   L      L 

exists  between  the  parameters  which  describe  the  incident  wave.  For  waves 
of  finite  amplitude  the  relationship  is  more  complicated  but  nevertheless 
exists. 

A  dimensional  analysis  of  the  physical  paramaters  indicated  in 
Equation  (1)  gives  the  following  dimensionless  parameters. 


e!zJH£„f|2»«>>l>H_r  Ji )  (3) 


pga2£     pga2£       L   a  2a  p  /gha^ 

As  pointed  out  by  Garrison  and  Chow  [5]  the  last  term  on  the  right  hand 
side  of  Equation  (3)  represents  the  ratio  of  the  Froude  number  to 
Reynolds  number  and,  accordingly,  is  an  indicator  of  the  ratio  of  the 
viscous  to  gravity  forces.  If  this  number  is  small,  corresponding  to 
large  scale  flows  with  small  viscosity,  it  may  be  assumed  that  this 
parameter  may  not  be  important  since  the  flow  would  be  controlled  primarily 
by  gravitational  forces.  In  this  case  Equation  (3)  may  be  written  as 


x  max    y  max  _  f  /27ra  h  H  \  ,.v 

pga2£     pga2£     I  L   a  2a  I 
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However,  as  noted  previously,  the  wave  period  may  be  used  in  place  of  the 
wave  length,  in  which  case  Equation  (4)  would  appear  in  alternate  form  as 

x  max    y  max  _  f  /gT2  h  H  \  .~. 

pga2£    pga2£     I  h   a  2a  I 

Equation  (5)  may  also  be  written  in  a  third  form  as 


Fx  max  or  Fymax  =  f  I  gT^  gT^  h  \  (6) 


pga2£     pga2£       h    H 


In  the  experijTental  program  an  attempt  is  made  at  determining  the  corre- 
lation of  the  dimensionless  parameters,  on  the  right  hand  side  of 
Equations  (4),  (5)  and  (6),  with  the  force  coefficients. 

The  complete  analytical  solution  to  the  wave  interaction  problem  is 
quite  complex.  However,  by  use  of  Morrison's  Equation  [4],  an  approxi- 
mate representation  may  be  carried  out.  For  this  purpose,  the  velocity 
potential  associated  with  the  incident  wave  is  written  as 

#  =  H  g  Cosh  k  (yfh)   ^^  _  at)  (?) 

2  a  Cosh  kh 

where  k  =  2tt/L  denotes  the  wave  number  and  a  =  2tt/T  denotes  the  frequency. 
From  Equation  (7)  the  horizontal  component  of  velocity  and  acceleration 
is  obtained,  respectively,  by  differentiation  as, 

u  m   _  H  gk  Cosh  k(y+h)  sin(kx  _  at)  (8) 

2  a  Cosh  kh 

.=Hgk  Cosh  k  (yfh)  ^^.^  (9) 

2    Cosh  kh 
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Ihe  horizontal  force  acting  on  the  cylinder  is  written,  according  to 
Morrison's  Equation,  as 

C, 

Fx  =  —  p  2a£  u2  +    (1  +  Cm)7ra2Jtu  (10) 

where  C.  denotes  the  drag  coefficient  and  C     denotes  the  added  mass 

d  J  m 

coefficient.  However,  for  cases  where  the  particle  motion  is  small  in 
comparison  to  the  cylinder  diameter  separation  effects  are  small  and  it 
is  possible  to  disregard  the  drag  contribution  to  the  total  force  and 
write  Equation  (10)  as 

-*-  -  (1  •  +  C  )  *k  5-  Cos  (  at)  (n) 

pga2£        m     2  Cosh  kh 

Using  the  definition  of  k  and  a,  Equation  (11)  may  be  written  in  terms 
of  the  dimensionless  parameters  indicated  in  Equation  (4)  as 


^^ =  (1  +  c)  2^2 ^^         (12) 

pga2MH/2a)        m   L  Cosh  (^  §■) 

ij  a 


The  vertical  component  of  the  force  could  be  expressed  in  a  form 
similar  to  Equation  (10) ,  but  the  vertical  component  of  acceleration  on 
the  bottom  is  zero.  It  is  supposed,  therefore,  that  the  vertical 
component  of  force  is  associated  with  the  lift  force  only.  Accordingly 
the  vertical  force  may  be  written  as 


F  =  £■  u2  C^  2a£  (13) 


2 
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in  which  CL  denotes  the  lift  coefficient  and  u  the  wave  induced  particle 
velocity  on  the  bottom.  Substituting  for  the  velocity  from  Equation  (8) , 
Equation  (13)  may  be  written  as 


Fy 2  n       j2Tia 

2, 


pgaz£ 


Sinfot_         (14) 


2a    Sinh  2kh 


Equations  (12)  and  (14)  contain  two  coefficients,  namely  the  added  mass 
coefficient/  C  ,  and  the  lift  coefficient,  CT .  Values  for  those  coeffi- 

m  Li 

cients  which  account  for  both  the  bottom  and  free  surface  are  not 
available  in  the  literature.  However,  if  the  water  is  relatively  deep 
the  free  surface  effect  should  be  small.  In  this  case  the  exact  potential 
flow  value  for  a  circular  cylinder  in  contact  with  the  rigid  bottom  has 
been  calculated  by  Garrison  [8]  as 

C  =  2.29. 

m 

For  the  same  geometry,  Daltman  and  Helfinstein  [9]  have  calculated  the 
lift  coefficient  by  an  approximate  method.  Their  value  is 


°L 


=  4.. 


These  values  of  the  two  coefficients  may  be  used  in  Equations  (12)  and 
(14)  to  obtain  approximate  expressions  for  the  horizontal  and  vertical 
force  coefficients.  It  should  be  noted,  however,  that  the  assumptions 
underlying  these  approximate  relationships  are  met  only  when  the  scale 
of  the  waves ,  i  „e . ,  the  wave  length  and  water  depth  are  large  in  compari- 
son to  the  cylinder  diameter.  Moreover,  since  linear  wave  theory  has  been 
used  to  determine  the  velocity  and  acceleration,  the  validity  is  also 
restricted  to  small  amplitude  waves. 
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III.   DESCRIPTION  OF  APPARATUS  AND  EXPERIMENTAL  PROCEDURE 

A.  APPARATUS 

In  order  to  measure  the  wave  forces  acting  on  the  horizontal  cylinder, 
experiments  were  carried  out  in  a  small  wave  channel.  The  basic  wave 
channel,  as  shown  in  Figure  2,   consisted  of  three  sections;  a  deep  water 
section,  a  shallow  water  section  and  a  transition  section,  each  fifteen 
inches  wide  with  a  rectangular  cross  section.  The  overall  channel  length 
was  sixty-four  feet. 

Three  quarter  inch  exterior  plywood  sheets  were  used  for  the  channel 
sides.  Vertical  rigidity  was  provided  by  bracing  the  plywood  sheets  with 
2  x  4 '  s  at  various  intervals ,  every  four  feet  in  the  shallow  water  section , 
every  foot  in  the  deep  water  section  and  every  two  feet  in  transition 
section.  Two  2  x  4's  were  attached  to  the  top  of  the  vertical  studs  and 
extended  the  length  of  the  channel. 

Since  the  bottom  of  the  deep  water  section  was  resting  on  a  concrete 
foundation,  only  a  single  section  of  plywood  was  used  for  the  channel 
floor.  A  double  floor  was  used  in  the  shallow  and  transition  sections. 
This  was  done  to  provide  extra  strength.  The  double  floor  was  constructed 
from  two  plywood  sheets  and  two  and  a  half  inch  separators.  The  sides  of 
these  sections  were  bolted  together,  through  the  spacing  of  the  double 
bottom,  with  three-eighths  threaded  rod.  Steel  angle  spacers  were  placed 
across  the  channel  top  to  maintain  the  width  dimension. 

Prior  to  assembly  all  wooden  sections  exposed  to  water  were  water- 
proofed. The  channel  interior  was  then  given  two  coats  of  Morewear 
Vitri-Glaz  1320A.  Dow  Corning  Sealant  780  was  applied  to  all  seams  and 
joints. 
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1.  Deep  Water  Section 

Several  methods  of  increasing  wave  height  and  steepness  are 
available.  One  of  the  most  common  methods  is  the  generation  of  waves  in 
a  deep  water  regime  and  allow  them  to  transition,  by  shoaling,  to  a 
shallow  water  regime.  By  the  proper  choice  of  depths  and  shoal  slope,  the 
transition  can  be  such  that  there  is  no  reflection  from  the  shoal,  and  the 
resultant  shallow  water  waves  are  in  a  near  breaking  condition.  It  was 
this  method  of  increasing  wave  height  that  was  chosen.  As  shown  in 
Figure  2,  the  deep  water  section  was  sixteen  feet  long  and  four  feet  deep. 
A  wave  generator  was  located  two  and  three-quarters  feet  from  the  closed 
end  of  the  channel.  In  order  to  generate  regular  waves  in  the  channel  a 
paddle  type  wave  generator  was  used.  An  aluminum  plate,  adequately 
reinforced,  was  hinged  at  the  wave  channel  floor  and  attached  by  a  pin 
connection  at  the  top  to  a  driving  rod,  Figure  3.  One-eighth  inch  teflon 
sheeting,  four  inches  by  fifty-one  inches  was  located  along  the  vertical 
sides  of  the  plate  to  prevent  leakage  from  one  side  to  the  other.  In  the 
space  behind  the  wave  paddle  the  wave  motion  was  damped  by  use  of  baffle 
plates.  A  two  horsepower  variable  speed  drive,  with  an  output  speed 
range  of  twenty  to  one  hundred  and  eighty  revolutions  per  minute,  was 
mounted  on  the  top  of  the  wave  channel  and  fitted  with  a  six  and  three- 
quarter  inch  radius  face  plate.  The  face  plate  was  provided  with  ways 
so  that  the  wave  generator  driving  rod  could  be  adjusted  to  any 
eccentricity  between  zero  and  eight  inches,  Figure  4. 

2.  Transition  Section 

As  previously  noted,  the  wave  height  was  increased  by  shoaling 
deep  water  waves.  To  accomplish  this  a  14.14  foot  long  ramp  was 
constructed  and  installed  between  the  deep  and  shallow  sections.  This 
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FIGURE  3.   WAVE  GENERATOR  PADDLE 
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provided  for  the  transition  from  deep  water  to  shallow  water.  In  order 
to  reduce  the  amount  of  reflection  from  the  transition  ramp,  the  slope 
was  set  at  1:6.  The  change  in  depth  due  to  the  shoal  was  two  feet. 
3.  Shallow  Water  Section 

The  shallcw  water  section  had  a  depth  of  two  feet  and  a  length 
of  forty-two  feet.  The  test  module  was  located  twenty  four  feet  from  the 
transition;  thus,  allowing  enough  length  for  the  waves  to  reach  a  fully 
developed  state  following  the  transition,  before  reaching  the  test  section. 
At  the  end  of  the  channel  was  located  a  variable  slope  beach  for  purposes 
of  dissipating  the  wave  motion.  The  beach  consisted  of  metal  shavings 
held  between  two  sheets  of  perforated  metal.  These  were  separated  by 
two  inch  wooden  spacers.  A  solid  sheet  of  one-eight  inch  aluminum  was 
attached  to  the  bottom  of  the  beach  by  three  and  five-eighths  inch 
separators.  On  the  exposed  surface  of  the  beach  was  attached,  parallel 
to  the  wave  fronts,  one  by  one  inch  aluminum  angles.  The  maximum  slope 
of  the  beach  was  1:7. 

a.  Test  Module 

In  order  to  provide  a  means  of  locating  the  test  cylinder  and 
measuring  the  wave  forces  acting  on  it,  a  separate  test  module  was 
designed.  The  test  module,  consisting  of  the  circular  cylinder  along  with 
its  own  walls  and  floor,  was  constructed  as  an  integral  system  in  order 
that  it  may  be  removed  from  the  channel.  This  concept  allowed  for  the 
maintenance  of  close  tolerances  during  manufacture  and  provided  for  easy 
access.  To  accept  the  test  module,  the  wave  channel  sides  and  floor  were 
recessed  and  made  of  plexiglass. 
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The  four  inch  diameter  plexiglass  cylinder  was  suspended 
between  the  walls  of  the  test  module  by  use  of  cantilever  beams  and 
adjusted  to  approximately  one-sixteenth  inch  above  the  floor,  Figure  5. 
It  was  known  a.  ptvloKi   that  this  sixteenth  of  an  inch  clearance  would 
cause  a  high  velocity  jet  under  the  cylinder  due  to  wave  action.  To 
prevent  this,  a  flexible  plastic  barrier  was  installed  in  the  gap.  The 
barrier  was  held  in  place  by  '0'  ring  material  pressed  into  0.007  inch 
slot  in  the  cylinder  and  the  test  module  floor. 

In  order  to  suspend  the  cylinder  a  small  distance  off  the 
wave  channel  floor  and  at  the  same  time  measure  forces  exerted  on  it  by 
surface  waves,  cantilever  beam  mounts  were  used  as  supports  as  shown  in 
Figures  6  and  7.  Bulkheads  were  fixed  in  the  cylinder  at  approximately 
two  inches  from  each  end,  and  the  fixed  ends  of  the  cantilever  beams 
bolted  on  these  bulkheads.  The  free  ends  of  the  two  beams  were  allowed 
to  protrude  slightly  beyond  the  end  of  the  cylinder  and  into  the  plexi- 
glass test  module  walls.  The  free  ends  were  supported  by  small  self- 
aligning  ball  bearings  pressed  into  the  plexiglass.  Both  beams  were 
fitted  with  strain  gages  and  waterproofed  using  BIH  Barrier  '  C f , 
Figure  8.  One  of  the  beams  with  its  largest  cross  sectional  dimension 
horizontal  was  used  to  measure  the  vertical  force  and  the  other  with  its 
largest  cross  sectional  dimension  vertical  was  used  to  measure  the 
horizontal  force.  Calibration  tests  shewed  the  cross  coupling  to  be 
negligible. 

This  design  presented  some  unique  problems.  The  cantilever 
beams  had  to  be  flexible  enough  to  allcw  measurement  of  the  forces  and 
stiff  enough  so  that  the  natural  frequency  was  large  in  comparison  to 
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FIGURE  6.  HORIZONTAL  FORCE  CANTILEVER 
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FIGURE  7 .  VERTICAL  FORCE  CANTILEVER 
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FIGURE  8.      FORCE  CANTILEVER 
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the  excitation  frequency.  The  characteristics  of  the  strain  gage  bridge 
and  amplif  ier  determined  the  minimum  strain  necessary  to  produce  a 
measurable  output.  The  strain  resulting  on  a  cantilever  from  a  given 
load  is 

6W. 


where 


**i 


e  =  strain 

P  =  load 

L  =  cantilever  beam  length 

E  =  modulus  of  elasticity 

b  =  cantilever  beam  base  dimension 

h  =  cantilever  beam  height  dimension 


(15) 


The  relationship  for  the  spring  constant  of  a  cantilever  beam  is 


Ebh-? 
K  =  —  (16) 


By  substituting  Equation  (15)  into  (16)  a  relationship  for  the  spring 
constant,  in  terms  of  the  strain  and  cantilever  dimension  is  determined 

ka   p3 
K  =  ZLl (17) 

e3E2b2h£ 

It  is  apparent  that  the  maximum  stiffness  and,  consequently,  the  maximum 
natural  frequency  is  obtained  by  making  b,  h  and  e  as  small  as  possible. 
The  minimum  value  of  b  was  set  by  the  width  necessary  to  mount  the  strain 
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gages.  The  minimum  value  of  e  was  determined  by  the  maximum  sensitivity 
of  the  amplifier/recorder.  The  maximum  stiffness  was  then  achieved  by 
making  h  as  small  as  possible.  However/  Equation  (15)  shows  that  the 
beam  length  decreases  with  b  and  h  as  expressed  by 

eEbhJ 
^   6P 

Thus,  the  lower  limit  on  h  was  finally  determined  by  setting  a  reasonable 
itujiimum  value  for  the  beam  length. 

In  order  to  prevent  the  cylinder  from  rotating,  an  offset  arm 
was  connected  from  the  horizontal  force  cantilever  to  the  test  module. 
A  Farber  bearing,  AMS  1K7,  was  used  as  a  wheel  at  the  test  module  end  of 
the  arm  to  reduce  friction. 
4.  Wave  Height  Probe 

There  are  two  types  of  probes  in  common  usage  for  measuring  wave 
height;  these  are  commonly  known  as  a  capacitance  and  a  resistance  type 
probe.  The  resistance  type  utilizes  two  vertical  wires  placed  a  small 
distance  apart;  the  amount  of  wetting  of  the  wires  varies  the  resistance 
between  them.  This  active  element  or  "variable  resistance"  is  wired 
into  one  leg  of  a  Wheatstone  bridge  and  the  signal  read  out  on  an 
amplifier/recorder.  Experience  has  shown,  however,  that  the  resistance 
type  probe  is  extremely  sensitive  and  well  suited  to  measuring  small 
amplitude  waves  but  is  rather  nonlinear  at  larger  amplitudes. 

The  capacitance  type  wave  height  probe,  although  less  sensitive, 
tends  to  be  linear  over  large  wave  amplitudes  and  was,  therefore,  used 
as  the  primary  method  of  sensing  wave  height  in  the  present  study.  The 
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probe  consists  essentially  of  a  single  insulated  wire  which  acts  as  a 
variable  capacitor;  the  wire  acts  as  one  plate  and  the  water  the  second 
plate.  As  the  water  surface  rises  and  falls,  the  capacitance  varies 
accordingly. 

The  probe  used  in  the  present  study  is  shown  in  Figure  9.  Located 
at  the  bottom  of  the  foil-like  support  member  was  a  three-eighths  inch 
diameter  acrylic  rod.  Hie  sensing  wire  was  terminated  at  the  bottom  in 
the  acrylic  rod  and  at  the  top  at  a  nylon  block.  The  tip  of  the  acrylic 
rod  was  made  removable  in  order  to  facilitate  changing  the  sensing  wire. 

A  number  of  wire  types  were  considered  and  tested  prior  to 
selecting  number  30  A.W.G.  wire  with  polythermaleze  insulation.  This 
wire  provided  the  capability  of  measuring  wave  heights  to  less  than  one 
thirty-second  of  an  inch  with  good  accuracy. 

The  electronic  circuit  utilized  with  this  probe  was  obtained 
in  schematic  form  from  Hewlett  Packard  and  is  shown  in  Figure  10 .  A 
Hewlett  Packard  carrier  amplifier  350  was  used  with  the  probe  circuit. 

B.  TEST  PROCEDURE 

With  the  channel  filled  to  the  desired  level,  the  test  section  and 
wave  height  probe  were  set  in  place.  In  each  instance  the  probe  was 
immersed  to  a  depth  of  nine  inches.  The  amplifiers  were  then  adjusted 
to  balance  the  bridge. 

With  the  use  of  a  transverse  mechanism  the  calibration  of  the  wave 
probe  was  accomplished.  The  probe  immersion  was  varied  over  a  seven 
inch  range  in  half -inch  increments.  The  gage  output  was  measured  both 
on  the  Hewlett  Packard  recorder  and  one  channel  of  a  Brush  recorder. 
This  provided  the  initial  calibration  information.  Prior  to  and  after 
each  set  of  runs  the  calibration  was  checked  at  several  points. 
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The  force  cantilevers  were  calibrated  by  exerting  a  load  on  the 
cylinder  by  means  of  a  series  of  weights  and  a  wire  pully  arrangement  as 
shown  in  Figure  11.  The  force  was  transmitted  to  the  cylinder  through  a 
series  of  tapped  holes  located  around  the  circumference  at  three,  nine 
and  twelve  o'clock.  Calibration  information  was  obtained  from  the  Bush 
recorder  output  by  loading  the  cylinder  in  half-pound  increments.  As 
with  the  wave  height  probe,  calibration  was  checked  prior  to  and  after 
each  series  of  runs. 

With  the  calibration  completed,  the  wave  generator  speed  was  adjusted 
for  the  desired  wave  length.  The  wave  height  probe  was  located  an  integer 
number  of  wave  lengths  frcm  the  cylinder.  Its  output  was  connected  to 
the  horizontal  grid  of  an  oscilloscope.  A  second  probe  was  then  located 
directly  over  the  cylinder  and  its  output  connected  to  the  vertical 
oscilloscope  grid.  By  observing  the  scope  pattern  the  first  probe  could 
be  adjusted  to  be  in  phase  with  the  waves  at  the  cylinder.  The  second 
probe  was  then  removed  and  a  data  run  commenced.  During  each  run  the 
wave  generator  drive  eccentricity  was  varied  from  a  half-inch  to  eight 
inches,  or  to  a  point  where  the  waves  broke  in  the  channel. 
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FIGURE  11.  CANTILEVER  CALIBRATION  ARRANGEMENT 
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IV.   PRESENTATION  OF  RESULTS  AND  CONCLUSIONS 

A.   EXPERIMENTAL  RESULTS 

As  discussed  in  Theoretical  Analysis,  non-dimensional  wave  forces  can 
be  represented  as  functions  of  relative  water  depth  (h/a) ,  relative  wave 
length  (2-na/L)   or  period  parameter  (h/gT2)  and  relative  height  parameter 
(H/2a) .  In  order  to  ascertain  the  effect  of  varying  these  parameters  on 
the  horizontal  and  vertical  components  of  wave  force  a  test  program  was 
established  which,  for  a  given  run,  held  water  depth  and  period  constant 
while  varying  the  wave  height.  Four  series  of  runs  were  conducted  at 
different  water  depths  corresponding  to  the  relative  water  depths  (h/a) 
of  9.0,  7.0,  5.5  and  4.0.  For  each  of  these  water  depths,  the  relative 
wave  length  (2i\a/L)   was  varied  over  a  range  from  0.06  to  0.48. 

Due  to  the  difficulty  in  presenting  all  the  wave  force  traces,  a 
series  of  representative  traces  is  shown  in  Figures  12,  13,  14  and  15. 
These  traces  represent  three  wave  lengths  from  the  three  deeper  water 
depths  (h/a  =  9.0,  7.0  and  5.5)  and  two  wave  lengths  for  the  shallowest 
water  depth  (h/a  =  4.0) .  For  each  of  the  eleven  runs,  four  wave  heights 
arranged  in  increasing  order  are  shown.  The  uppermost  trace  represents 
wave  height,  the  middle  vertical  force  and  the  bottom  trace  is  horizontal 
force. 

The  wave  height  traces  show  that  the  waves  were  generally  sinusoidal 
in  nature.  This  was  especially  true  at  the  shorter  wave  lengths.  At  the 
greater  wave  lengths  the  wave  trace  departed  frcm  the  sinusoidal  form  as 
the  wave  height  was  increased.  As  can  be  seen  in  Figure  12,  for  the 
greater  lengths,  as  the  height  was  increased  the  wave  was  characterized  by 
sharp  peaks  and  long  flat  troughs. 
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Before  discussing  the  wave  force  traces  a  brief  discussion  of  the 
mechanism  involved  in  the  wave  forces  is  in  order.  In  Bernoulli's 
equation  two  contributions  to  the  pressure  acting  on  the  surface  of  the 
cylinder  are  recognized,  an  unsteady  contribution  associated  with  the  time 
derivative  of  the  potential  to  the  first  power  (vertical  forces)  and  a 
term  associated  with  the  fluid  velocity  squared.  The  velocity  squared 
term  tends  to  be  symmetric  about  the  y  axis  and,  consequently,  contributes 
little  to  the  horizontal  force.  The  inertia  term  is  therefore  dominant, 
in  the  horizontal  forces,  and  increases  linearly  with  wave  height  H/2a. 
Examination  of  the  representative  traces  shows  that  the  horizontal  force 
varied  nearly  sinusoidally  with  wave  height  for  all  cases  except  when  the 
wave  length  became  very  large.  In  the  case  of  very  long  waves,  the  hori- 
zontal force  peaks  are  separated  by  a  null  region  which  corresponded  to 
the  wave  trough.  The  maximum  value  of  the  horizontal  force  coincided  with 
the  maximum  water  particle  acceleration  associated  with  the  incident  wave, 
Equation  (9) . 

For  the  vertical  forces,  there  are  two  regimes,  one  in  which  the 
inertial  forces  dominate  and  a  second  in  which  the  velocity  squared  term 
dominates.  The  inertial  contribution  to  the  pressure  produces  both  posi- 
tive and  negative  values  of  vertical  force.  This  term  tends  to  be  linear 
in  wave  height  H/2a  and  have  a  maximum  value  of  uplift  when  the  cylinder 
is  under  the  wave  trough.  The  velocity  squared  term  produces  only  upward 
forces  due  to  the  non-symmetrical  flow  about  the  cylinder.  This  component 
of  vertical  force  is  a  function  of  the  wave  height  squared  (H/2a) 2  and 
peaks  occur  at  the  points  of  maximum  wave  induced  velocity  occurring  at 
both  crest  and  trough  of  the  wave.  Examination  of  the  vertical  force 
traces  shows  that  for  the  longer  waves  the  force  was  regular  and  had  a 
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frequency  twice  that  of  the  incident  wave.  The  double  frequency  was 
caused  by  the  velocity  inaxinrums  at  the  crests  and  troughs.  The  maximum 
occurring  under  the  wave  crest  in  the  longest  waves  was  apparently 
attributable  to  the  fact  that  the  velocity  under  the  crest  of  a  shallow 
water  wave  is  greater  than  at  the  trough.  Here  the  upward  force  was  lift 
dominated.  At  intermediate  wave  lengths  the  reverse  occurred.  As  the 
wave  length  was  shortened  the  inertial  effects  gradually  took  effect  and 
the  maximum  occurred  at  the  trough  where  the  two  components  were  additive. 
Finally,  at  very  short  wave  lengths  the  inertial  force  became  dominant 
and  the  traces  showed  a  sinusoidal  behavior;  the  lift  force  associated  with 
the  velocity  squared  term  in  Bernoulli's  equation  became  negligible. 

1.  Horizontal  Forces 

As  noted  previously,  the  horizontal  force  was  nearly  sinusoidal  for 
the  entire  series  of  data  runs.  It  is  therefore  possible  to  characterize 
the  force  by  the  maximum  value.  In  order  to  examine  the  variation  in  the 
horizontal  force  with  varying  wave  height,  plots  were  constructed  of  the 
non-dimensional  force  coefficient  fx  (fx  =  F    .Vpga2£)  versus  the 
relative  wave  height  H/2a.  These  wave  force  versus  wave  height  plots  are 
presented  in  Appendix  B.  It  is  apparent  from  the  plots  that  horizontal 
force  coefficient  is  characterized  by  a  linear  variation  with  wave  height. 
Inasmuch  as  the  variation  of  the  force  coefficient  was  linear  with  H/2a, 
for  a  given  wave  length,  it  is  possible  to  characterize  the  horizontal 
force  by  the  dimensionless  parameter  fx/ (H/2a) ,  the  slope  of  the  wave 
force-height  plots. 

The  parameter,  F  _/[pga2 (H/2a)]  is  presented  as  a  function  of  the 
relative  wave  length  (2iTa/L)  in  Figure  16.  This  figure  also  presents  the 
result  predicted  by  Morrison's  analysis,  Equation  (12)  .  For  the  larger 
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values  of  relative  wave  length  (2-na/L)   the  experimental  results  and 
theory  are  in  good  agreement.  This  is  generally  as  expected  because  for 
the  larger  values  of  2-na/L  the  waves  are  well  represented  by  linear  theory 
and  the  effect  of  the  free  surface  is  reduced  due  to  the  greater  water 
depth.  It  may  be  recalled  that  these  were  the  assumptions  upon  which  the 
Morrison  equation  is  based.  However,  for  smaller  values  of  2-na/L,   the 
experimental  results  are  consistently  high.  This  is  not  unexpected  as  the 
waves  have  become  shallow  water  waves  and  are  not  well  represented  by  the 
linear  wave  theory.  Also,  in  the  case  of  the  smaller  depths  the  effect  of 
the  free  surface  becomes  important. 

2.  Vertical  Forces 

As  previously  discussed,  the  vertical  force  may  display  two  regimes, 
one  inertia  dominated  and  one  dominated  by  the  force  associated  with  the 
velocity  squared  term  in  Bernoulli's  equation.  The  assumption  having  been 
made  in  Equation  (13)  that  the  inertia  forces  would  be  negligible,  the 
vertical  forces  would  be  a  function  of  the  velocity  squared  only  and  hence 
a  function  of  (H/2a) 2 .  It  is  difficult  to  characterize  the  force  vari- 
ations with  wave  height  by  a  single  term  as  two  types  of  flow  exist. 
However,  having  made  the  assumption  that  the  force  was  a  function  of 
velocity  squared,  it  was  decided  to  use  fy/(H/2a)2  as  the  characteristic 
quantity. 

Values  of  F  ^/pga2£(H/2a)  2  were  plotted  against  the  relative  wave 
length  (2iTa/L)  and  presented  in  Figure  17.  The  values  of  this  coefficient 
predicted  by  a  lift  force  only  are  also  shown  in  this  figure.  The  agree- 
ment between  the  theory  and  experimental  results  appears  to  be  good  for 
the  greater  water  depth  (h/a  =  9.0)  .  However,  as  the  water  depth  was 
decreased  the  agreement  became  poorer  and  the  experimental  data  displayed 
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higher  force  values  than  predicted  by  Equation  (14) ;  apparently  owing  to 
the  effect  of  the  free  surface  and  the  non-linear  effects  of  the  incident 
wave.  Also,  as  the  wave  length  decreased,  or  the  parameter  2Tra/L  increased, 
it  may  be  noted  that  the  agreement  is  poor.  This  is  apparently  due  to  the 
fact  that  the  inertia  component  becomes  dominant  and  this  contribution  to 
the  force  was  neglected  in  Equation  (13) . 

As  waves  are  often  characterized  by  period  rather  than  wave  length, 
Figure  18  presents  the  force  coefficient  as  functions  of  a  relative  period 

(VgT2) . 

B.   CONCLUSIONS 

From  the  experimental  results  the  following  conclusions  are  considered 
warranted: 

1.  The  horizontal  forces  on  the  cylinder,  resulting  from  the  incident 
gravity  waves,  increased  linearly  with  wave  height. 

2.  The  vertical  force  showed  two  regimes,  one  inertial  dominated  and 
one  velocity  squared  dominated.  The  inertial  force  regime  occurred  at 
shorter  wave  lengths,  while  the  velocity  squared  term  predominated  at  the 
longer  wave  lengths. 

3.  The  maximum  vertical  force  was  always  smaller  than  the  maximum 
horizontal  force  and  the  downward  vertical  force  was  nearly  zero. 

4.  For  greater  water  depths  the  theoretical  predictions  appeared  to 
give  valid  results  for  both  horizontal  and  vertical  forces  in  the  short 
wave  length  range. 

5.  Horizontal  forces  became  independent  of  the  parameter,  2-na/L,   for 
greater  wave  lengths. 
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